Inter-Landau level magnetoexcitons in bilayer graphene 



(N 

o 

(N 

> 
O 

(N 



o3 



T3' 

c: 
o. 

>: 

00, 

try 

p; 

ov 
o: 

(N 



X 



Judit Sari 1 and Csaba Toke 2 
1 Institute of Physics, University of Pecs, H-7624 Pecs, Hungary and 
2 BME-MTA Exotic Quantum Phases "Lendulet" Research Group, 
Budapest Univ. of Technology and Economics, Institute of Physics, Budafoki ut 8., H-llll Budapest, Hungary 

(Dated: November 26, 2012) 

If bilayer graphene is placed in a high perpendicular magnetic field, several quantum Hall plateaus 
are observed at low enough temperatures. Of these, the a xy = 4ne 2 /h sequence (n / 0) is explained 
by standard Landau quantization, while the other integer plateaus arise due to interactions. The low- 
energy excitations in both cases are magnetoexcitons, whose dispersion relation depends on single- 
and many-body effects in a complicated manner. Analyzing the magnetoexciton modes in bilayer 
graphene, we find that the mixing of different Landau level transitions not only renormalizes them, 
but essentially changes their spectra and orbital character at finite wave length. These predictions 
can be probed in inelastic light scattering experiments. 

PACS numbers: 71.35.Ji, 71.70.Di, 71.70.Gm 



I. INTRODUCTION 

Bilayer graphene^ two coupled hexagonal lattices of 
carbon atoms in the Bernal stacking^ of graphite, is 
a two-dimensional zero-gap semiconductor with chiral 
charge carriers with Berry's phase 2tt, having a roughly 
parabolic dispersion at low energies about the corners of 
the hexagonal first Brillouin zone4 These facts are testi- 
fied by its unusual integer quantum Hall effect^ 5 - (IQHE) , 
featuring a double step in the ladder of the Hall conduc- 
tance in a strong perpendicular magnetic field B±, first 
observed by Novoselov et alX This double step, 8e 2 /h in- 
stead of the common 4e 2 /h for spin and valley degenerate 
Landau levels, is due to the degeneracy of the n = 0, 1 
Landau orbitals, 6 - - — The gap at the integer quantum Hall 
effect at filling factor v = phc/eB± = ±4, ±8 has been 
recently measured with great accuracy^ and the excita- 
tions of the IQHE states in the long wavelength limit 
have also been observed by infrared absorption^ and 
Raman spectrocopyAi Further broken symmetry states 
have been observed^— in the central Landau band at 
v = 0, ±1, ±2 and ±3, and by careful tilted-field measure- 
ments it has been shown that they arise predominantly 
from many-body effects, i.e., from quantum Hall ferro- 
magnetism (QHF). 16 Quantum Hall states with broken 
symmetry have also been found in the n = — 2 Landau 
levelii, and there is also some evidence for a fractional 
quantum Hall plateau. 

The eightfold degeneracy of the central group of Lan- 
dau levels is at best approximate, because the Zeeman 
energy is unavoidably present. While the latter is rather 
small on the characteristic scale of the interaction en- 
ergy, a perpendicular electric field can be applied to bias 
the two layers; 8 ' 14 ' — which causes an energy difference 
between the two valleys. The competition of the on-site 
energy difference between the layers and interactions may 
result in interesting physics, especially at v = 0i 14 ' 31 ~ — 

If the chemical potential is in the gap between Lan- 
dau bands, the low-energy excitations are bound particle- 
hole pairsj^Zr— called magnetoexcitons .As the net charge 



of such an excitation is zero, taking appropriate linear 
combinations one obtains eigenstates of the total mo- 
mentum. In such states the hole and the particle are 
bound by the attractive Coulomb interaction, thus form 
a dipole with a separation of q£ B at center-of-mass wave- 
vector q, where Ib = \J h/eB is the magnetic length. 
These modes determine the transport gap in the q — > oo 
limit. Some of these modes couple to circularly polar- 
ized light while others may be observable in Raman 
scattering experiments ! 11 ' 42 ' 43 

For monolayer graphene Yang et al-^ studied the 
intra-Landau level excitations of quantum Hall ferromag- 
netic states, and Iyengar et a/4 5 - and Lozovik et al^ dis- 
cussed the inter-Landau level excitations in detail both 
for the IQHE and QHF states. Infrared absorption data 
by Jiang et and Henriksen et al£& indicate a contri- 
bution from many-body effects to the inter-Landau level 
transitions in these states. 

In bilayer graphene Henriksen et alJ£ have found that 
fitting single-body parameters does not fully explain the 
observed cyclotron resonance; Deacon et alM> and Zouet 
al£^ have found a significant particle-hole asymmetry, 
whose origin is still debated^ The wave-vector depen- 
dence of the excitations have not been observed so far. 
Theoretically, the intra-Landau level excitations in bi- 
layer graphene are well understood. At odd integers in 
the central Landau band (y — ±1,±3) the degeneracy 
of the n = 0, 1 Landau orbitals causes fluctuations with 
an in-plane electric dipole character, which gives rise to 
unusual collective modes . 16 ' 52 ' 53 At even integers in the 
central Landau band (v = ±2, 0) orbital degeneracy does 
not play a similar role, and the intra-level excitations are 
still magnetoexcitons i 34 ' 54 To complement these studies, 
here we address the issue of the inter-Landau level exci- 
tations of bilayer graphene in the quantum Hall regime. 

Our paper is organized as follows. In Sec.|H]we review 
the tight-binding model of bilayer graphene, and the ba- 
sic facts concerning its Landau levels and orbitals. Our 
goal is to systematically explore the range of applicabil- 
ity of subsequent simplified models, which neglect several 
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parameters of the Slonczewski-Weiss-McClure 55 (SWM) 
model, or account for them on the level of perturbation 
theory. We intend to add a few observations to the excel- 
lent studies available in the literature^ - — In Sec. IIIII we 
review the adaptation of the mean-field theory of magne- 
toexcitons to the case of bilayer graphene. In Sec. II VI we 
study the excitations of the IQHE states, and in Sec. |V] 
those of the QHF states. We conclude in Sec. I VII with 
an outlook on experimental connections. 



II. LANDAU LEVELS AND ORBITALS 

Each layer of bilayer graphene consists of two sublat- 
tices, denoted A and B in the top layer and A and B 
in the bottom layer. In Bernal stacking 2 - two sublattices, 
A and B in our notation, are exactly above/below one 
another, while the A sites are above the center of the 
hexagons in the bottom layer, and B sites are below the 
centers of hexagons in the top layer. See Fig. [T] 




FIG. 1. (Color online) Bilayer graphene in Bernal stacking. 
The hopping parameters of the Slonczewski-Weiss-McClure 
model,— conventionally denoted 70, 71, 73 and 74, are also 
indicated. 

The low-energy physics of bilayer graphene can be ad- 
equately described by the tight-binding effective theo- 
ries that specialize the SWM models of graphite to the 
case of just two layers. In the vicinity of the valley cen- 
ters corresponding to the K (£ = 1) and K' (£ = —1) 
first Brillouin zone corners, this amounts to using the 
Hamiltonian£ 
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where n — p x +ip y and p = — iftV — eA, v — \/3ajo/2h 
10 6 m/s is the intra-layer velocity, v% — v / 3a73/2?i is the 
trigonal warping parameter, 71 is the inter-layer hopping 
amplitude, and V4 = V3a'j4 : /2h is a velocity parame- 
ter related to interlayer next-nearest neighbor hopping. 

= is the Zeeman energy (with g being the 

gyromagnetic factor and //g the Bohr magneton). This 



Hamiltonian acts in the basis of sublattice Bloch states 
[ipA> ^b'^a^b] in valley K and [^ig, ^U> V>B> ^] hi val- 
ley K 1 . Here 70 = jab = Iab ^ s the intra-layer hopping 
amplitude, 71 = is the interlayer hopping amplitude 
between sites above each other in the two layers. Further, 
73 = 7 A jj and 74 = 7^ = j B g are next-nearest neigh- 
bor interlayer hopping amplitudes, as shown in Fig.[TJ A' 
is the on-site energy difference between the dimer sites 
(A,B) and the non-dimer sites (B,A). Finally, u is the 
potential energy difference between the layers, which may 
arise, e.g., because of an applied perpendicular electric 
field E±. 

The Hamiltonian in Eq. (fT]) is block-diagonal in the 
valley index, which is conveniently described as a pseu- 
dospin. In the special case u — the system has SU(2) 
pseudospin rotation symmetry. In the theoretical limit 
Az —> this is raised to SU(4) symmetry. In this paper 
we will treat Az and u as small perturbations in com- 
parison to the interaction energy, i.e., we will work in the 
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< 1, 



< 1 limit. We set h = 1. 



For small momenta, p -C 71 /4w, the two low-energy 
bands of the Hamiltonian that touch each other at 
K and K' in the case of vanishing magnetic field can be 
attributed to a 2 x 2 effective Hamiltoniar>i£ 
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where m = 71 /2i> 2 , and H'^ acts on [ipA, ipg] in valley K 
and [i/>|>, "0a] in valley K' . The Landau levels (LL's) and 
Landau orbitals, respectively, of the two-band Hamilto- 
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Here n =^ —I is an integer, is an appropriate normal- 
ization factor, and r\ nq are the single-particle states in the 
conventional two-dimensional electron gas with quadratic 
dispersion in the Landau gauge A = yBx, 
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and H n is a Hermite-polynomial. 

The n = 0, 1 orbitals are degenerate in the u —> 
limit, and they have a layer polarization for £ = ±1. 
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At realistic values of Az and u, the n — 0,1, £ = ±1, 
a =t, | states form a quasidegenerate band we will call 



the central Landau level octet. Notice that A 



(n) 



(n)| 



l/\/2 for ii — y 0. 



and |S 

We would like to determine how neglecting 73, 74, and 
A' changes the single-body orbitals of is Eq. ([1}, 
and how much these differ from the simplified two-band 
model H'^ in Eq. ([2]) . The values of the SWM parameters 
for bilayer graphene was estimated by a combination of 
infrared response analysis and theoretical techniques by 
Zhang et alM- They found 



— = 0.133, — 

7o 7o 



A' 

0.1, — 

7o 



0.006. 



(7) 



These ratios are based on 70 = 3.0 eV. While somewhat 
greater values of 70 are also available in the literature^ 
we use these values for a robustness analysis. For the 
particle-hole symmetry breaking term we use 



74 
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0.063 



(8) 



from the recent electron and hole mass measurement by 
Zou et al.r^ which is slightly greater than the value in 
Ref . [13. 

With 7r = -jf^a, [a, cl*] = 1, the Hamiltonian can be 
expressed in terms of these Landau level ladder opera- 
tors. Then the eigenstates of can be calculated nu- 
merically. Fig. [5] shows the overlap of the Landau orbitals 
with the "ideal" limit 73 = 74 = A' = as the SWM pa- 
rameters 73,74, A' are tuned from zero to their literary 
values [Eqs. (J7) to ©] for the central in = 0, 1) and 
the two pairs of nearby (n — ±2, ±3) Landau levels. At 
small magnetic field B = 0.1 T trigonal warping alone 
significantly changes the orbitals from their ideal limit. 
Switching on 74 hardly affects the central levels, but for 
n > 2 it changes the electron and hole pairs (+n, — n) 
differently, as expected from this electron-hole symmetry 
breaking term. Finally, the inclusion of A' hardly affects 
the orbitals. These changes, however, are already small 
at modest fields (B = 1 T), and are further suppressed 
as the experimentally relevant range (B « 10 T) is ap- 
proached. Thus neglecting the 73, 74, A' SWM param- 
eters is justified in the high magnetic field range where 
quantum Hall experiments are typically performed. 

As the two-band model in Eq. @ applies for small mo- 
menta, and the low-index Landau orbitals have a small 
amplitude at high momenta, the two-band model is ex- 
pected to be valid for the lowest few Landau levels. The 
Landau orbitals of the two-band model have a large over- 
lap with those of the four-band model in the "ideal" 
limit 73 = 74 = A' = 0: 1, 0.9995, 0.9992, 0.9987 for 
n = 0,1, ±2, ±3, respectively. We conclude that using 
the Landau states of two-band model in Eq. ([2]) instead 
of these of the four-band model in the lowest-energy 
Landau bands does not introduce further inaccuracy be- 
yond the neglect of 73, 74 and A'. Therefore, we take H'^ 
as our starting point. 
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FIG. 2. (Color online) The overlap of the Landau orbitals 
with the "ideal" limit, 73 = 74 = A' = 0, as the SWM 
parameters 73,74, A' are gradually tuned from zero to their 
literary values [Eqs. to ©] for the lowest-energy Landau 
levels. For the effect of 73 in the two-band model see Ref. IBTj . 



III. MAGNETOEXCITONS 

When the Fermi energy is in a Landau gap (e.g. at 
filling factor v = . . . , -12, -8, -4, 4, 8, 12, . . . in bilayer 
graphene), the integer quantum Hall effect 4 occurs in 
samples with moderate disorder! 1 - Quantum Hall states 
also occur at other integer filling factors because the 
exchange interaction favors symmetry-breaking ground 
states called quantum Hall ferromagnets; single-body 
terms such as the Zeeman energy play a secondary role. 
QHF's emerge at odd integer fillings in two-component 
systems even if the Zeeman energy is tuned to zero. 
This observation straightforwardly generalizes for SU(n) 
systems^ 

Because of the clear separation of the filled and empty 
Landau bands in the (mean-field) ground-state, the ex- 
citations of both classes of quantum Hall systems are 
described in the same way. The relevant low-energy exci- 
tations are magnetoexcitons, which are obtained by pro- 
moting an electron from a filled Landau band to an empty 
bandi^Zr— These neutral excitations have a well-defined 
center-of-mass momentum Q. They approach widely sep- 
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arated particle-hole pairs in the Q oo limit. The lat- 
ter limit determines the transport gap unless skyrmions 
form.'"' 6 Magnetoexcitons are created from the ground- 
state by operators 3 ^— 

p 

where N — (n,£, a) (N' — (n',£',er')) specifies the Lan- 
dau band where the particle (hole) is created and A is 
the area of the sample. 

Magnetoexcitons carry spin and pseudospin (valley) 
quantum numbers, as derived from the particle and hole 
Landau bands involved. While the projections S Z ,P Z of 
the spin and the pseudospin are always good quantum 
numbers, their magnitudes 5* and P are well-defined only 
for ground-states that are spin or pseudospin singlets, re- 
spectively. 

It is common practic o 37 i 45 to define the quantity 

l z = \n\ - \ri\, (10) 

and consider it the "angular momentum quantum num- 
ber" of the exciton. We emphasize that l z is exactly 
conserved by the electron-electron interaction only in the 
Q — > limit, where it is related to angular momentum. 
The emergence of this quantity is best seen in the two- 
body problem of the negatively charged electron and the 
positively charged hole, as discussed in the Appendix. At 
any finite wave vector transitions with different l z may 
mix. 

In the low magnetoexciton density limit the interaction 
between magnetoexcitons is neglected. The mean-field 
(Hartree-Fock) Hamiltonian of magnetoexcitons is well 
known from the literature ^2r— and so is its adaptation 
to spinorial orbitalsi 34 i 45 

H (NN')'y L i) - 
(G|%^(Q)y*Lv,(Q)|0) - 5 N x6 N ,fi,{0\V\0) = 

= ^NN^N'N' (Ena(t7 ~ ^n'a'i'cr' + A(n, 7l')) + 

+ <S(Q) + <S(Q), (ii) 



where N = (n,£,cr), etc., and Snn> = 5 aa >5^i5 nn ' ■ The 
first term of the r.h.s. is the single-body energy difference 
of the N and N' states, which includes the wave vector 
independent exchange self-energy difference of the two 
states. While the exchange self-energy itself is infinite 
for any orbital, its difference between two states, 

X N > N = J ?2^2^'^(P)' ( 12 ) 

N'M — Xmn) , (13) 

M filled 

is finite. (We will define I^^(p) soon.) This is a pecu- 
liarity of the two-band model of bilayer graphene, which 
is replaced by an appropriate renormalization procedure 
for monolayer graphene 4 ^, and for the four-band model 
of bilayer graphene. 34 

The next term is the direct dynamical interaction be- 
tween the electron and the hole: 

<S(Q) - -/ ^e*<«*«>/£* (q). (14) 

This term is diagonal both in spin and pseudospin, oc 
Sd-o-S^Sd-io-'Sf;,^ , but not in Landau orbital indices. Fi- 
nally, the last term in Eq. (fTTj) is the exchange interaction 
between the electron and the hole, 

<^(Q) = 2^r^S'(Q). (is) 

which is oc S a(T iS^iSaa'S^, , thus couples transitions that 
conserve the spin a and the valley £ of the electron and 
the hole individually. Sometimes we will call it the RPA 
contribution, as it is related to particle-hole annihilation 

and recreation processes. Notice R[nn'](Q) vanishes in 
the Q — > limit. 

We have used the notation 



'S(P) = V(q) K?^4-M^ 1 M^^ rjV2||Ar , | (p)F |JVl||Ar , | (p)+ 



i(™2) Am) An{) Ana) 



R (™ 2 ) r("i) nW) r(™2) j-i* / \ 771 / \ 

B e 2 % B (,' 1 % ^1^1-2,1^1-2, (P)- f, |iVi|-2,|7V;|-2(P) 



|(™ 2 ) p(ll) R ("l) Am) 



(™ 2 ) Anx) A (n[) R (n 2 ) , 
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V2 



e -? 2 4/4 if n > > nj else F NN ,((i) = F* N , N (-q). 



Here N — 2 = (n — 2, £cr), | N | = (\n\, f, c), and -FVjv(q) is related to the Fourier transform of rj nq (r) in Eq. ([6]), 
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and L™(z) is an associated Laguerre polynomial. The 
difference between the intralayer Coulomb interaction 
V s (q) = 2ne 2 / (Aireoe r q) and the intralayer one, V D (q) = 
e~ qd V (q) where d ~ 0.335 nm is the distance between 
the layers, is neglected as a first approximation. The 
and numbers correspond to the spinorial struc- 
ture of the single-body states in Eq. A 1 = A = 1, 

B 1 = B° = and A^ l) = sgn(n)S^ n) = l/y/2 for n > 2. 
Notice that 



E, 



(NN') 



:q) 



i?P;j(Q) = R{S,](Q) = (-l)"+"' +ft+fi '^^(Q), 

(17) 



E 



(NN') 
(NN 1 ) * 



Q) = ^ ) ) (Q), (16) 



?(N'N), 
L (N'N) y 



which follow from the similar properties of I^^?{p). 

The q — > limit of the magnetoexciton dispersion de- 
termines the many-body contribution to the cyclotron 
resonance. In the case of the ordinary two-dimensional 
electron gas, traditionally implemented in semiconduc- 
tor heterostructures and quantum wells, this mush van- 
ish by Kohn's theorem. 60 This theorem, however, relies 
upon the simple quadratic form of the kinetic energy op- 
erator, thus it is not applicable in monolayer and bilayer 
graphene. Thus it is possible that the sum of the ex- 
change self-energy constant Xjvjv an d the diagonal part 
of the dynamical interaction i?^^,j(Q) does not vanish 
in the Q — > limit. 

The mean field Hamiltonian matrix H^^,} (Q) in gen- 
eral mixes transitions among different electron-hole pairs, 
restricted only by conservation laws. Landau level mix- 
ing effectively screens the interaction. Sometimes the 
magnetoexciton spectra are obtained using a screened 
model interaction instead of the bare Coulomb, not let- 
ting LL transitions mixj 51 i 59 We believe such an approach 
is suitable in the q = limit, where an additional quan- 
tum number l z also restricts LL mixing, and for intra-LL 
modes. At finite wave vector the mean-field theory with 
LL mixing removes spurious level crossings in the excita- 
tion spectra and provides insight into the orbital struc- 
ture of the excitations. Technically, however, the infinite 

H(nn')(Q) ma t r i x needs to be truncated. 



IV. INTEGER QUANTUM HALL STATES 

We first consider the states where the chemical poten- 
tial is between two orbital Landau bands. This occurs at 
filling factor v = . . . , -12, -8, -4, 4, 8, 12, . . . in bilayer 
graphene. Together with S z and P z , the magnitude of 
the spin S and of the pseudospin P are quantum num- 
bers. (In the Az — > 0, u — > limit an SU(4) classification 
is also possible.) 

With the hole (n!) and the electron (n) Landau or- 
bitals fixed, the sixteen possible transitions belong to four 
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FIG. 3. Magnetoexciton modes at the v — —4 integer quan- 
tum Hall effect in bilayer graphene. Only the optically rele- 
vant spin- and pseudospin-conserving modes are shown. In 
the q —¥ limit, which is probed by purely optical ex- 
periments, the transitions in the infinite sequence at fixed 
lz — \n\ — jn'| may mix; at q > 0, the excitations result from 
the mixing of all l z sequences. The modes at v = +4 are 
obtained by particle-hole conjugation. 



classes: (i) a spin singlet, pseudospin singlet state: 

€(q) = ^EE*U(q). (is) 

(ii) A spin singlet, pseudospin triplet multiplet. The 
P z = member of this multiplet is 

(iii) A spin triplet, pseudospin singlet, which contains 
following the S z = state 

s ° 

(iv) A nine-member multiplet that is triplet is both spin 
and valley. Its S z — = P z member is 

^n'(q) ^EE^W^Uy^W- ( 21 ) 

Fig. [3] depicts these modes for v = —4. 
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The exchange interaction between the electron and the 
hole contributes only to states generated by {f2} nll /. 
In all other excitation modes the RPA term cancels (c.f. 
the sign alternation^ in modes ^]^" and ^J^) 

or is prohibited by quantum numbers. Thus in the ab- 
sence of accidental degeneracies, we expect a collection 
of nondegenerate excitations and one of fifteen-fold de- 
generate excitations; the latter is decomposed as 4, 7, 4 if 
A z > 0, u = and as 1, 2, 1, 2, 3, 2, 1, 2, 1 if A z , u>0. 

The extent interactions may mix transitions involving 
different Landau level pairs depends on the interaction- 
to-kinetic energy ratio, parametrized by 

P = — ~~~f /fr^c oc — 1 -=, (22) 

where e r is the relative dielectric constant of the environ- 
ment. LL mixing is eventually suppressed in the B — > oo 
limit just like for the conventional two-dimensional elec- 
tron gas. Realistically (10 T < B < 40 T, 1 < e r < 4), 
1 < /3 < 8; this is by no means a small perturbation. 

For q — the mean-field Hamiltonian mixes magne- 
toexcitons with different electron and hole Landau lev- 
els at fixed l z , and for q > it also mixes different l z 
subspaces; c.f. Fig. [3] Restricting LL mixing to a fixed 
l z subspace might give the impression that LL mixing 
is just a quantitative correction for the long wave length 
part of the lowest excitation curves, resulting in increased 
electron- hole binding energies . 46 ' 62 However, already the 
lowest excitations in the different l z sectors mix strongly 
at finite wave length. As the side panels of Fig. 0] show, 
the excitations have a large projection on the l z subspaces 
different from their own l z in the q — > limit, and may 
eventually be contained in one of the other subspaces for 
large q; this is an unavoidable consequence of the elim- 
ination of crossings by LL mixing. Thus in the rest of 
this paper we allow the mixing of transitions restricted 
by \l z \ < L and a maximum number M at each fixed 
l z . Fig. [4] shows transitions at v = 4 with L = 1 and 
M = 7, while Fig. [5] shows \v\ — 8 and 12. We will use 
this truncation in the spectra shown in the rest of this 
paper <£i 

The q — > limit of the magnetoexcitons are commonly 
probed by optical absorption and electronic Raman scat- 
tering. The selection rule o 41 ' 43 ensure that only the ^j^ , 
mode and the ^l^i , &w and &tii modes of the fifteen- 
fold degenerate curve are active, l z = ±1 is absorption 
and l z = in Raman. Particle-hole conjugation relates 
v = 4n to v = —An (n positive integer) with the sign of 
l z reversed. 



V. QUANTUM HALL FERROMAGNETIC 
STATES 

With an integer filling factor different from v = 
±4, ±8, ±12, . . . , a Landau band quartet (\v\ > 4) or 
octet < 4) is partially filled in the single elec- 



tron picture. The minimization of the interaction en- 
ergy results in gapped states which break either spin 
rotation or pseudospin (valley) rotation symmetry, or 
bothJ^^^^SS If either the Zeeman energy A z or 
the interlayer energy difference u is present, they deter- 
mine the sequence how the Landau levels are filled j 1 ^ The 
most convenient basis in pseudospin space may differ; we 
may introduce 



a n scrp = cos -fi n) £=i )<rp + sin -e l0 a ni5= _i iCTp , (23) 

9 
a n Acrp = sin -a n ^ =1 ^ p - cos -e l, *a ni5= _i i(T p. (24) 



With proper choice of 9 and <f>, these definitions include 
states of definite valley, bonding and antibonding states, 
or intervalley phase coherent states. Corresponding mag- 
netoexciton operators are defined in an obvious manner. 

In particular, if Az > u, the v = QHF state is 
ferromagnetic and the choice of the pseudospin basis is 
irrelevant (Fig.[S]). For v — ±2, both the n — and n = 1 
orbital Landau levels of identical spin and pseudospin 
are filled, where <f> and 9 are determined by electrostatics 
(Fig. [9}. For odd v, an interlayer phase coherent (0 < 9 < 
7r/2) state exists for sufficiently small enough u, which 
yields^ to a layer polarized state (9 = 0) at v = —3 
and v = 1, and to a sequence of states with partial or 
full orbital coherence^ 3 - at v = — 1 and v = 3. Notice 
that v — —3 and v = 3 are not related by particle-hole 
symmetry. This is best understood from Hund's rules^: 
at v = — 3 only one n — orbital band is occupied, while 
at v = 3 only one n = 1 orbital band is empty. At 
v = 3 Cote et at 53 showed that orbitally coherent states 
dominate the phase diagram, whose inter-LL excitations 
are beyond the scope of this study. The case of v = — 3 
is depicted in Fig. [11] 

Beyond S z and P z , the magnitudes P or S are quan- 
tum numbers at half-filling v = 0. All excitons include 
transitions between the non-central levels |n| > 2; the 
possibility of transitions from, to, or within the central 
Landau level octet depends on the ground state, which 
also resolves the transitions through the exchange self- 
energy differences; see Figs. [51 [51 and [TT1 

The excitations are grouped by their optical signature. 
Due to the small momenta of optical photons, valley flip- 
ping modes are optically inactive. In the q — > limit 
l z becomes a quantum number, and l z = ±1 applies for 
single-photon absorption^ 1 -, and l z = 0,±2 in electronic 
Raman processes 4 ^, with the l z = transitions being 
dominant. The angular momentum due to the helicity of 
the photons is transferre d 41 ' 43 entirely to the orbital de- 
gree of freedom. Optically inactive modes include Gold- 
stone modes associated with the broken symmetry (out- 
side the scope of our study) and generic dark modes. 




FIG. 4. (Color online) The excitations of the integer quantum Hall state at \v\ = 4 and B — 10T. The mixing of Landau 
levels is truncated at L — 1 and M = 7. Solid lines show the fifteen-fold degenerate excitations, which include three optically 
relevant S z — Pz — modes. Dashed lines show the spin and pseudospin singlets. Inset: spectra if Landau level mixing is 
neglected. Side panels: the weight of the definite l z projections in each curve in bottom-up order. The top row shows the spin 
and pseudospin singlets. 



It is known that two QHF ground states exist, 
a spin-polarized one and a valley (layer) polarized 
on o 14 ' 16 ' 31 ' 32 ' 34 — . Their respective range of validity is 
determined by the ratio of the Zeeman energy Az to the 
energy difference between the valleys, which in turn is 
related to the potential difference u, (In fact, layer and 
pseudospin can de identified within the central Landau 
level octet.) For concreteness, we are discussing the fer- 
romagnetic state. Here the magnitude of the pseudospin 
P is a good quantum number of the excitations. See 
Fig. [5] for the transitions that span the Hilbert space of 
the mean-field Hamiltonian. 

With the electron (n) and hole (n') Landau levels fixed, 
the S z = — 1 transitions consist of a pseudospin triplet^ 5 - 

*nA|.,Ti'St' ^nSi,n'Ap T/f (^nSi,n'St ~ ^nAUn'Af) ' anc ^ 

a singlet 65 (^ nSit s t + ^ nA i,At) ■ This group con- 
tains the intralevel transitions among the n = 0, 1 Lan- 
dau bands; the Goldstone modes associated with the spin 
rotational symmetry breaking should be in this subspace. 
However, our approach is not appropriate for the descrip- 
tion of Goldstone modes even for even filling factors, as 
we will discuss below. 

The S z = 1 triplet and singlet of the 



reversed 



component 



structure, 



t 

nSf.n'Ap 



j_ (V 

V2 v " 



nSt,n'Si 



t 

nA1,n'Al I ' 



t 

nAt,n'Sl' 

and 



respectively, contains 



inter-LL transitions only. 

The S z — sector consists of (i) two triplets, 



\F/^ \T/^ J_ I \[/t _ \T/t A 

an ^ ^liAi,n'Sp *nSJ.,n'A4.! 7/f (^nSi,n> 



- \J>' 

Si nAl,n'Al 

the RPA terms does not contribute to, and (ii) 



and 



two singlets, -J- (*^ sti „, st -r- - nAt , nMt 

^(*nS4,n'S4. + *Lu,n'Al.)' which &K mixed by 
the RPA term. Careful inspection reveals, how- 
ever, that the two pseudospin singlets (ii) always 
appear in the mean-field Hamiltonian on equal foot- 
ing, e.g., the ^ (*2St,iSt + *2At,iAt) transition is 
indistinguishable on the mean-field level from the 

by 



2, Si + ^ lAi-2,Ai 



transition. This follows 



p(2,i) _ p (i -2) 

^(2,1) — ^(1,-2) ' 



(25) 

d(2,1) _ p(l -2) _ _d(1-2) _ _p(2,l) / 9R \ 
•"(2,1) — "(1,-2) — (2,1) — Xl (l,-2)' 

and the following easily provable identity of the exchange 
self-energy cost: 

A(n,n')+X n >o + X n >i-X ni0 -X ntl = A(-n',-n). 

(27) 

(For n = 1 or n' = 1 no sign change is necessary.) 
Eq. (|27|) simply expresses particle-hole symmetry, i.e., 
that the exchange self-energy cost of transitions related 
by particle-hole conjugation in a fixed component must 
be identical. See Fig. [7] for the transitions whose pom- 
parison yields Eq. (|27l) . 

The RPA terms are the same in each diagonal and off- 
diagonal position among equivalent transitions, thus they 



FIG. 5. (Color online) The excitations of the integer quantum 
Hall state at \v\ = 8 and \v\ = 12, at B = 10T. 



select the even and the odd linear combinations in group 
(ii). The even combination, ^™,(q) defined in Eq. ([TBI) , 
gets an RPA enhancement, while the RPA cancels from 
the alternating sign combinations, making it energeti- 
cally equivalent to the P z = element of the triplets (i) . 
Thus, eventually, the S z = 0,P Z =0 sector contains a 
three-fold denegerate curve and a nondegenerate mode. 

Each of the four multiplets in the S z = sector con- 
tains a P z — mode, which is active in electronic Ra- 
man or IR absorption. Here the mixing of Landau levels 
results in more widely separated modes. The optically 
active excitations are shown in Fig. [5] with LL mixing 
taken into account. Notice that the l z = 1 and l z = — 1 
transitions have an equal weight in all modes, consistent 
with the particle-hole symmetry at v — 0. 

In the S z = — 1 sector we find spin waves. Neglect- 
ing Landau level mixing, they give rise to a gapless and a 
gapped intra-LL modes^ and a sequence of higher inter- 
LL modes; each of these is raised by the Zeeman energy 
and split by the valley energy difference in turn. The in- 
teraction, however, mixes these excitations, thus a clear- 
cut classification into intra-LL and inter-LL is no longer 
possible. Level repulsion unavoidably lowers the formerly 
gapless modes. This effect yields apparently negative ex- 
citation energies at small wavelength. Goldstone's theo- 
rem, however, ensures that a gapless spin-wave mode is 
associated with the breaking of the spin rotation sym- 



FIG. 6. Magnetoexciton modes at the v = quantum Hall 
ferromagnetic state in bilayer graphene. Only the optically 
relevant spin- and pseudospin-conserving modes are shows. 
In the q — s> limit, which is probed by photon absorption or 
electronic Raman, the transitions in the infinite sequence at 
fixed l„ — \n\ — \n'\ may mix; at q > 0, the excitations result 
from the mixing of all l z sequences. 

metry. Consequently, the seemingly negative energies of 
the lowermost curves with a large intra-LL component 
is an artifact of the combination of Hartree-Fock mean- 
field theory and LL mixing. The same anomaly occurs 
for monolayer graphen o 45 i 46 , but it is less apparent when 
the particle-hole binding energy is plotted. 

B. v = ±2 

The v = ±2 state breaks the spin and pseudospin ro- 
tational symmetries as the ground state fills the n = 
and n = 1 orbitals of the most favorable spin-pseudospin 
component, S f. We restrict the discussion to spin and 
pseudospin preserving excitations. See Fig. [9] for the pos- 
sible transitions. 

It is easy to check that the mean-field Hamiltonian ma- 
trix is identical to the one at v = 0. For the — n — > n 
transitions [n > 2 integer) this holds because the occu- 
pancy of the central Landau level octet is irrelevant as 

X-nfl + X- n i — X n fi — X n i = 0. (28) 

The octet of — [n + 1) — > n and — n — > (n + 1) transi- 
tions gives rise to two quartets of equivalent transitions 
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FIG. 7. Imagine two transitions: (I) a particle in level n' 
is promoted to level n of a valley-spin component (left) of 
which the n = 1 Landau band is filled, and (II) — n — > —n 
of a component (right) of which the n = 1 band is empty. 
Particle-hole conjugation of (I) is the promotion of a hole from 
Landau level — n' to — n, i.e., (II); the excitation energies must 
be the same. The consequence is Eq. (|27|l ; the formal proof 
is straightforward. 



by Eqs. (|25H27|) . While at v = the S t and A f tran- 
sitions of the former group bundle with the S 1 and A 1 
transitions of the second group, now the S f transition 
of the first group bundle with the A "f, S 1 and A 4- tran- 
sitions of the second group. The spectrum is still the one 
in Fig. |SJa). The orbital projection of the modes differs, 
c.f. Fig. QJJfor v = —2. At v = +2 the sign of l z changes 
in all projections, which determines the helicity of the 
absorbed and inelastically scattered photons. 



C. f = -3 

At v = —3 the only occupied band is (0,5* f). The 
n = orbitals are selected by their small energetic 
advantaged The states in Eq. (|2"3")l progress from the 
layer balanced limit 9 = tt at u = to the layer po- 
larized state 6 = 0; this limit is achieved about u = 
u c « 0.001e 2 /e^s, which is only 0.082 meV at B = 20 T. 
For < u < u c there is interlayer phase coherence*^ 
Thus magnetoexcitons exist on both side of u c ; the 
amount electrostatics raises energy of the pseudospin- 
flipping modes w.r.t. the pseudospin conserving modes 
saturates at u = u c . Both the spin and the pseudospin 
symmetries are broken resulting in a three-fold Goldstone 
modej£2 Further, Barlas et ah— showed that at finite u 
there is an instability to a stripe ordered phase with a 
rather small critical temperature. Our analysis below 
applies only below this temperature. See Fig. [TT] for the 
transitions that span the Hilbert space of the mean-field 
Hamiltonian. 

Because of the degeneracy of n = 0, 1 orbitals, states in 
central Landau level octet at odd integer fillings involve 
fluctuations with in-plane electric dipole character J£ The 
consequent collective modes have been studied in detail 
by Barlas et al2^ and Cote et al^ As we do not han- 
dle such dipolar interactions, we have omitted the pre- 
dominantly intra-LL lowest curve from the spectra, and 



we have checked that the inter-LL excitation modes we 
keep contain the 0—^1 magnetoexcitons with a negli- 
gible weight. Reassuringly, we always got a weight less 
than 0.1%. 

As the n = orbital is filled with S t electrons in the 
mean-field ground state, for fixed electron (n) and hole 
(n r ) Landau levels the exchange self-energy cost of the 
S t transition is higher than those of the other compo- 
nents. Also, in the S f component an intralevel — >• 1 
transition is possible, which mixes with higher S f tran- 
sitions; see Fig. [TT] for the restrictions on the possible 
transitions at this filling. The other three components, 
on the other hand, occur symmetrically in the mean-field 
Hamiltonian; c.f. Fig. [TT] One can change basis from the 
excitons of type S 1, A f and A 1 to 



i 



j, J 2 



(29) 



2* 



(30) 



\j>t r — _L f\T/t i v|/t I flyt A 



(31) 



The RPA contribution cancels from ^ nn i(q) and 
^(q), which give rise to doubly degenerate excita- 
tions, ^^/(q) has an RPA contribution. Its mixture 
with the distinguished S f excitations produces singlet 
excitations. In higher energy excitations, on the other 
hand, the weight of the — > 1 transition of the S f com- 
ponent becomes extremely small, thus the equivalence 
of components will be approximately restored, yielding 
threefold quasi-degenerate and nondegenerate curves. 

See Fig. [TH for the dispersion of the active excitations. 
The small graphs show that the projection to definite l z 
subspaces changes abruptly at nonzero wave vector. No- 
tice that at higher energies > le 2 /(47reo£r) the double 
degenerate curves occur in the vicinity of a nondegen- 
erate one, indicating the approximate restoration of the 
equivalence of components in this limit. 



D. 



Just like at \v\ = 3, the state at v = — 1 is not the 
particle-hole conjugate of the state v = 1. We do not 
study v = — 1 because of the relavance of orbitally coher- 
ent states^ 2 - 

At v = 1 the interlayer coherent and layer polar- 
ized QHF states both have magnetoexcitons. The pos- 
sible transitions in the mean-field Hamiltonian are ob- 
tained trivially by raising the Fermi level by four levels 
in Fig. [TT] Now the S i electrons are distinguished by 
the possibility of an intralevel — > 1 transition, and their 
higher self-energy. The other three components occur 
symmetrically in the mean-field Hamiltonian. The argu- 
ment is similar to the case of v = —3; for the two 1—^2 
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FIG. 8. (Color online) (a) Excitation spectra of the quantum Hall ferromagnet at v — 0, ±2 at B — 10 T. Only the optically 
relevant S z = P z = modes are included, (b) The weight of the shown excitations on the definite l z subspaces for v = in 
bottom-up order. For degenerate curves the weight is summed. The quantum numbers of the q — > limit are indicated. 
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FIG. 9. The spin- and pseudospin-conserving magnetoexciton 
modes at the v = — 2 quantum Hall ferromagnetic state in 
bilayer graphene. The modes at v = +2 are obtained by 
particle-hole conjugation. 



transitions of spin-f electrons and the —2^1 transi- 
tions of the spin- \ electrons one uses -Rn'oj = — ^(i'o) 2 ^ 
and Eq. (|2T|) . When subsequent l z = ±1 transitions are 
included, the —n — > (n + 1) transition of A \, is equivalent 
to the —(n+1) —> n transition of S f and A f, while the 





-,1 ;-l 



12 34 5 12 34 50 12 34 5 12 3 



FIG. 10. (Color online) The projection of the excitations at 
v = —2 on the definite l z subspaces. For degenerate curves 
the weight is summed. The spectrum at v = ±2 is identical to 
the one in the left panel of Fig. O the letters refer to the same 
curves. For v = +2 the sign of the +/„ and — l z projections 
are interchanged w.r.t. v — —2. 



— n — > (n + 1) transition of A t and S f is equivalent to 
the — (n+1) — > n transition of A The convenient basis 
change is similar to Eqs. (l29l) to (|3Tj) . with S f replaced 
by A\.. 

There is only a slight difference between the inter- 
Landau level excitation spectra at v = — 3 and v = 1: 
the exchange self-energy cost of the distinguished transi- 
tion (S t at v = — 3 and S 1 at v = 1) relative to the three 
equivalent ones is X u iq — X n Q at v = —3 and — X„i 
at v = 1. As this difference is already small for the low- 
est transitions and then decreases, we omit the v = 1 
spectrum; its difference from Fig. I12f a) is comparable to 
the line width. 
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FIG. 11. The optically spin- and pseudospin-conserving mag- 
netoexciton modes at the v = —3 quantum Hall ferromagnetic 
state in bilayer graphene. Notice that the S \., A~\ and A J. 
transitions occur symmetrically, i.e., the same electron and 
hole Landau levels are allowed and the same self-energy cost 
is picked up from the exchange with the filled levels. 



VI. CONCLUSION 

We have calculated the inter-Landau level magnetoex- 
citons in the integer quantum Hall states as wells as 
the quantum Hall ferromagnets at filling factor v — 
—3, ±2, 0, 1 of bilayer graphene. We have found that the 
spinorial structure of the orbitals together with the en- 
hanced electron-hole exchange interaction effects in this 
multicomponent system gives rise to rather complex dis- 
persions; these are related both to the shape of the 
Fourier transform of Landau orbitals and to the elimina- 
tion of crossings by the mixing of Landau levels, which 
is significant because the scale of the Coulomb energy is 
comparable to the cyclotron energies. 

The q — > limit of the excitation can be 
probed by optical absorption 4 ^ and electronic Ra- 
man experiments i 1 1 1 43 Unlike for the conventional two- 
dimensional electron gas with quadratic dispersion^ the 
excitations in a quantizing magnetic field do acquire an 
interaction shift; the magnitude of such a shift is one of 
our experimental predictions. The wave vector depen- 
dence of the magnetoexcitons can be probed in inelastic 
light scattering experiments, 42 where the nonzero mo- 
mentum is taken away, or is provided by, an impurity or 



a phonon. The mixing of transitions that involve differ- 
ent Landau levels is strong in the experimentally acces- 
sible range. This mixing smoothens the dispersion rela- 
tions via level repulsion, and at finite wavelength causes 
a strong mixing of the modes that have distinct angular 
momenta in the zero wavevector limit. In particular, we 
have found that the anticrossings due to Landau level 
mixing result in a number of undulations in the magne- 
toexciton dispersions, whose van Hove singularities must 
give a strong signal. The elaboration of the relevant Ra- 
man scattering processes is delegated to further studies. 

On the theory side, we have found that the classifica- 
tion of magnetoexciton modes at finite wave vectors by 
the "angular momentum quantum number" is problem- 
atic, that the screening effects — which we have handled 
via Landau level mixing — is significant, and that this 
framework is not quite suitable for studying the Gold- 
stone modes such as spin waves in the symmetry-breaking 
quantum Hall ferromagnetic states. 
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Appendix A: The two-body problem in bilayer 
graphene 



Let us consider the Hamiltonian 



Ho = — 



1 



2m \tt 







(A 
o 



t\2 



1 

2m 







-u{t x -r 2 ) 



■Wi,yi Pi = -iWi 





(Al) 

(e/c)A(ri) belongs 



ifrS/2 + (e/c)A(r 2 ) belongs to 
2 /er. (We have fixed the valley 



where n t = p hX 
to the electron and p 2 
the hole, and u(r) = e 
of both the electron and the hole. A valley- independent 
interaction is assumed because any deviation from this is 
small in the ratio of the lattice constant to the magnetic 
length.) Introducing center-of-mass and relative coordi- 
nates (X,Y) = ri ^ r2 , (x, y) = ri — r 2 and momenta 

(P x ,Py) = ri + r 2 , ipx,P v ) = P1 2 P2 , and separating the 
center-of-mass motion by the canonical transformation 
U = e lXy , we obtain 



H' 2 



( Cl {C\f ^ 

(Cl) 2 (cl) 2 

c\ c 2 

V c\ {c ] _f / 



-u(r— zxP), 



(A2) 



where the independent harmonic oscillators C± are de- 
fined as 



^4 



(p x - zf ) ± ijpy - i|) 

2V2 



(A3) 
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FIG. 12. (Color online) (a) Excitation spectrum of the quantum Hall ferromagnet in the v = —3 at B — 10 T. Only the 
optically relevant S z = P z ~ modes are included, (b) The weight of the shown excitations on the definite l z subspaces in 
bottom-up order. For degenerate curves the weight is summed. 



In complete analogy to the case of the monolayer—, the 
eigenstates of the kinetic part of H' 2 are 



I S{nJ)<t>\ n+ \ t \ n _\_2 \ 

<t>\n+\,\n-\ 

S{n + )S(n_)(t)\ n+l _2 y \ n _\-2 
V S(n + )0|„ + |_2,| n _| J 



(A4) 



in terms of two-dimensional harmonic oscillator eigen- 
states, 



</w(r) 



2TTlr 



(A5) 



Above S(n) = is if n — or n — 1, else it is sgn(n). 
Thus the kinetic energy operator clearly commutes with 
the operator 



L = C\C+ -ClC_ 



/-2 0\ 





\ 2/ 



which returns l z as an eigenvalue. This operator, how- 
ever, commutes with the complete H' 2 only in the P — » 
limit. Thus we can regard the electron-hole bound state 
as a two-dimensional harmonic oscillator with clockwise 
(C-|-) and anti-clockwise (C_) excitations placed in an ex- 
ternal confinement potential. A nonzero center-of-mass 
motion breaks the rotational symmetry of this confine- 
ment and starts to couple the . n _ states. 

Here we closely follow Iyengar et al£^ The two-body 
Hamiltonian in Eq. (|A1|) can be written is terms of center- 



of-mass and relative coordinates and momenta as 



2^ T + 2p * 



(Py ~ 



2{p y - X) 



2m 



2m 



(2P xPx - 2{ Py - X)(Py - X)) 
(P X (Py - x) + 2 Px (Py - X)) 



+ ^(Px(Py-X) + (p y -X)P x + 

+Px(Py ~x) + (P y - x)p x ) - u(r), 



(A6) 



where m x — \(<J x ® 1 ± 1 ® cr x ), and = — \{<J y <E> 1 ± 
1 <S> (J y ). With the application of the canonical transfor- 
mation U = e lXy , we obtain 



fr, m t (V 2 , , 2 x 2 „ , : 



2m y 

^ (yx - 4p' x p' y ) 
2m J 

m~ 



(A7) 



{-yPy - PyV + Px X + X Px) 

u(r — z x P'). 



Eq. [jA"2)) follows by substituting Eq. (|A"3|) into Eq. (|AT| . 
Notice that ^ in Eq. (|A2[) is independent of the trans- 
formed center-of-mass coordinates. Thus P' is conserved; 
in original variables this corresponds to P — er x B. 
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